ABNORMAL PRESSURE DETERMINATION USING 
NUCLEAR MAGNETIC RESONANCE LOGGING 

FIELD OF THE INVENTION 

[0001] This invention relates generally to abnormal formation pressure 
determination in a borehole, and more specifically to such abnormality determination 
using a method and apparatus for making pulsed nuclear magnetic resonance (NMR) 
measurements of earth formations. Preferred embodiments are directed to abnormal 
formation pressure determination via NMR-derived porosity, and logging-while- 
drilling (LWD) or measuring-while-drilling (MWD) applications. This patent 
application claims priority of provisional patent application Ser. No. 60/459,490, filed 
April 1, 2003, which is incorporated herein by reference. 

BACKGROUND OF THE INVENTION 

[0002] Pore and fi-acture pressures are the major factors affecting the success of 
drilling operations. If pressure is not properly evaluated, it can lead to drilling 
problems such as lost circulation, blowouts, stuck pipe, hole instability, and excessive 
costs. Unfortunately, formation pressures can be very difficult to quantify precisely 
where unusual, or abnormal, pressures exist. Abnormal pressure is a subsurface 
condition in which the pore pressure of a geologic formation exceeds, or less 
fi-equently is below, the expected hydrostatic formation pressure. Abnormally high 
formation pressures typically occur when impermeable rocks such as shales are 
compacted rapidly during their deposition, so that their pore fluids cannot always 
escape and must then support the total overlying rock column. Such pressures may 
occur as shallow as a few hundred meters below the surface or at depths exceeding 
thousands of meters and can be present in shale/sand sequences and/or carbonate 
sections. 

[0003] The occurrence of overpressured zones has classically been determined 
well after the overpressured zone has been drilled, sometimes with disastrous results. 
Abnormal pressure can cause a well to blow out or become uncontrollable during 
drilling, resulting in overwhelming increases in drilling/completion costs. Therefore, 
the earliest possible determination of the presence of abnormally overpressured 
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subsurface(s) is a paramount piece of information that is of acute interest to the 
drilling engineer and the operating company. In particular, knowledge of formation 
pressure during the drilling enables the well operator to make preparations such as 
increasing the weight of the drill mud column in order to maintain well control and to 
prevent a blow out (or to resist fluid loss), as zones of overpressure (or underpressure) 
are penetrated by the drill bit. 

[0004] Generally, as an overpressured formation is approached, there are marked 
differences in degree of compaction and porosity of various formation layers. Also, 
differences in the minerals composition of interstitial fluids typically occur, and the 
formation pressure may begin to rise and approach or even exceed the bottom hole 
pressure, thus decreasing the bottom hole pressure differential. If properties, which 
are affected by such factors, are closely monitored and plotted with respect to depth, 
an abnormally pressured zone may be identified when a distinct deviation from an 
average line trending with depth is observed. To this end, the most commonly used 
methods for formation pressure prediction in logging-while-drilling (LWD) utilize the 
same phenomena: changes in shale porosity due to burial and deviation from "normal 
compaction" trend in abnormally pressurized zones. 

[0005] The most commonly used measurement techniques include resistivity and 
acoustic (sonic) logs that may be supplemented with gamma ray and/or SP 
(spontaneous potential) as shale indicators. In complex cases, a gamma-gamma 
density and neutron (pulse-neutron) logs are also used. Typically, such logging- 
while-drilling techniques are used to: (1) identify shale layers, (2) assess their 
porosity, and (3) account for possible complications, such as lithology changes due to 
diagenesis, water salinity variations, natural or drilling-induced fractures. The 
measured parameters, such as electric (resistivity), acoustic (travel time) and nuclear 
(gamma-gamma/neutron-gamma/pulsed neutron), are plotted versus depth and an 
established normal compaction trend is provided. Any deviation from this trend 
should be recognized on porosity logs and may be associated with a subsurface 
pressure anomaly. 

[0006] All conventional logs, however, suffer from one key weakness - estimated 
total porosity is not lithology-independent. As known, lithology-dependent methods 
give a correct answer in shales if: (1) the lithology is known; (2) there are no fracture 
components; and (3) the pore space fluid is pure water. However, uncertainties in 
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matrix composition, such as (1) mixed lithology and diagenetic alteration, (2) changes 
in pore water salinity (salt influence or fresh-water entrapment), (3) presence of 
hydrocarbons, and (4) natural/induced fractures cause calculation of formation 
pressure to become very challenging, if not impossible in complex reservoir cases. 
Most importantly, all known techniques have failed in cases of variations in water 
salinity (presence of salt dome) and fracturing (tectonic faulting/ drilling damage). 

[0007] The resistivity-based logging techniques are most effected by variations in 
matrix mineralogy. In normally pressured shales (hydro-pressure regime), there is a 
decrease of conductivity with depth of burial, due to compaction and associated 
porosity and water loss. This decreasing of shale porosity forms a "normal 
compaction line" trend in conductivity or resistivity versus depth plots. Generally, 
such a nomial compaction line can be established for a given area. Depending on 
various factors this line may shift right or left, but the slope remains essentially the 
same. Shale-reading points that fall to the left (higher conductivity area) of the 
normal compaction line are often associated with zones of abnormally high pressure. 
However, a basic assumption in the use of resistivity logs is that water salinity in 
shales remains the same, and therefore resistivity changes could be correlated to 
actual porosity variations. This assumption becomes invalid if there are variations in 
salt deposition, which may facilitate or impede water conductivity, thus rendering 
resistivity-based porosity values incorrect. 

[0008] The above problem is illustrated in Fig. 1, which is a resistivity log of an 
overpressured formation from Gulf-of-Mexico area. In the Fig. 1, the normal 
compaction trend is shown as soHd line, and resistivity data as black points. The 
deviation of resistivity at the depth of 2800m is associated with increased amount of 
clay-bound water (hence increased conductivity). In this well, two pressure 
compartment are identified: --2500-3200 m and 3400-3800 m. Reasons for this 
compartmentalization, however, cannot be determined from the resistivity log, since 
there is no information on the spectral distribution of the total porosity. Reasons for 
increasing resistivity at 3300 m are unclear also; they may be due to lithology 
changes, such as presence of dispersed carbonaceous material, diagenetic alterations 
in shales, or others. 

[0009] Acoustic logs are also effected by variations in matrix mineralogy. In 
normally pressured shales sonic logs generally show a decrease in interval transit time 
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with depth of burial. Since interval transit time is a function of porosity, this decrease 
indicates that the shale porosity decreases with depth. When depth is plotted on the 
linear y-axis versus shale interval transit time on the logarithmic x-axis, a straight line 
can be drawn trough the normal pressure points. Points that fall to the left of the 
normal line represent a zone of abnormally high formation pressure, as shown in Fig. 
2. Such acoustic log is less dependent on water conductivity, but strongly influenced 
by natural and drilling-induced fractures that may drive calculated porosity values too 
high. Moreover, both acoustic and density methods are strongly influenced by gas 
presence. 

[0010] Fig. 2 illustrates acoustic (sonic) log in an overpressured formation from 
Gulf-of-Mexico area. A normal compaction trend derived from the acoustic log is 
shown as the smooth line. Deviations to the left are due to increasing of shale 
porosity and hence are associated with undercompaction in overpressured formations. 
Top of pressure seal is observed at depth of about 3220 meters, where deviation of 
actual data points (circles) from the normal trend becomes apparent. Pressure seal is a 
gradual change in porosity related to pressure increase, also known as transitional 
zone. Overpressured formation zone marked by constant values in travel time, and 
started at approximately 4000 meters. 

[0011] Difficulties associated with the prior art are further related to the fact that 
often pressure predictions from various conventional logs disagree with each other. 
Fig. 3, for example, illustrates such disagreement in pore pressure estimates from 
resistivity and acoustic logs. Logging-while-drilling resistivity measurements 
(circles) show a long transitional zone in the 9,000-10,000 feet - interval, which is 
typical for damaged shaly pressure seal. Such a long transitional zone is probably 
associated with diagenetic alterations and corresponding variations in water salinity. 
In contrast, acoustic log (small dots in the midsection of the drawing) shows fairly 
steep changes, meaning that pressure seal is in good condition. Also, fracture gradient 
(dashed line), calculated based on pore pressure is in disagreement with leak-of test 
results (EMW - equivalent mud weight) line. Moreover, below 1 1000 feet 
calculations based on resistivity seriously overestimates formation pore pressure. 

[0012] Therefore, it is an object of the present invention to develop a method and 
system for real-time analysis of shales and accurate porosity-based formation pressure 
estimation to facilitate early warning of the existence of overpressure. In particular, it 
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is desirable for the new method to enable lithology-independent porosity 
measurements - that is, at least a portion of the extracted data is related only to shale 
micro-porosity, wherein diagenetic changes associated with shale composition are not 
interfering with the desired data. Additionally, it is desirable that at least a portion of 
the data can provide other information about a reservoir or shale zone. Additionally, 
it is desirable for the system to have mechanical strength and measurement sensitivity 
to withstand shock, vibration and erosion associated with drilling and to enable 
logging and measurement while drilling of the underground formation in most 
complex underground conditions. 

SUMMARY OF THE INVENTION 

[0013] To this end, in one aspect the invention is a nuclear magnetic resonance 
(NMR)-based method and system for real-time analysis of geologic formations 
including accurate porosity-based formation pressure estimation. In particular, the 
method and system of the present invention effectively resolve problems associated 
with prior art logging techniques discussed above, NMR's unique abiUty to 
discriminate porosity associated with clay-bound, capillary-bound and free fluid 
volumes allows the user to extract the signal only from shaly section and ignore the 
presence of "additional" porosity from induced fracturing or presence of silty 
sandstone material. Moreover, changes in pore water salinity are do not affect NMR 
measurements, and hydrocarbon signal can be recognized and accounted for using 
changes in the hydrogen index (HI) parameter. In general, since NMR logging used 
in accordance with the present invention is lithology-independent, it is not effected by 
matrix variations and provides accurate porosity estimation in mixed or even 
unknown matrix composition. 

[0014] In particular, it has been recognized by the present inventors that Ti 
(longitudinal polarization time) measurements are especially effective for while- 
drilling measurement of porosity. United States patent No. 6,051,973, the entire 
contents of which is incorporated herein by reference, describes a method of obtaining 
reservoir information based on Ti measurements. Unlike the T2 measurements that 
are sensitive to magnetic field changes (which may result, for example, from 
movement of the logging tool) and its associated dephasing effects, the Ti 
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polarization process is substantially immune to small field fluctuations characteristic 
to logging-while-drilling operations. In a preferred embodiment of the invention, 
these? motion-insensitive Ti measurements are used as a primary source of spectral 
porosity associated with clay-and capillary-bound water. Since an excessive porosity 
value characterizes shales in or overlying overpressured zones, these overpressured 
zones can be recognized using Ti measurements. 

[0015] Fundamentally, the proposed NMR approach utilizes the same principle as 
conventional porosity logs to provide advance identification of overpressure; namely, 
the overpressured subsurface shale intervals contain higher volumes of water within 
versus shale intervals following normal pressure trends. However, in accordance with 
a preferred embodiment, the method of the present invention provides an additional 
benefit of discriminating between total pore space volvmie, Free Fluid volume (FFI), 
Capillary Bound Water volume (BVI) and Clay Bound Water volume (CB W) of 
subsurface layers, while drilling the physical well. In accordance with one aspect of 
the invention, the direct continuous measurement of clay-bound water during the 
drilling process serves as definitive indication that an over-pressured zone has been 
entered. 

[0016] Accordingly, in a preferred embodiment, a method for determining 
abnormal pressure zones of a geologic formation is disclosed. In particular, in a 
preferred embodiment, a normal compaction trend is constructed using NMR-derived 
CBW content for non-consolidated formations or CBW and BVI for consolidated 
subsurface formations versus depth. Then, any deviations fi-om this normal 
compaction trend are determined using one of the industry-accepted techniques. Such 
deviations indicate a presence of shale sections with higher porosity that directly 
corresponds to an overpressured top seal or transitional zone. 

[0017] In accordance with different embodiments of the invention, various 
techniques may be used to determine pressure deviations fi^om the normal compaction 
trend. Such known techniques are described, for example, by (1) Hottmann et al. in 
"Estimation of Formation Pressures From Log-Derived Shale Properties," JPT 
{Journal Petroleum Technology) June, 1965; (2) Eaton in "The Equation for 
Geopressure Prediction From Well Logs," paper SPE 5544, 1975 Annual Technical 
Conference , Dallas Sept. 28 - October 1, 1975; (3) Rasmus et al. "A Framework To 
Estimate Pore Pressures in Real-Time," February 1990, Measurement- While-Drilling 



6 



NYJD: 1422032.1 



Symposium, Louisiana State University, Baton Rouge, LA; (4) Holbrook et al. in "A 
Petrophysical/Mechanical Mathematical Model for Real-Time Wellsite Pore- 
Pressure/Fracture-Gradient Prediction," Paper SPE 16666 presented at the 1987 SPE 
Annual Technical Conference, Dallas, Texas, September 1987; and (5) Bryant in "A 
Dual Shale Pore Pressure Detection Technique," paper SPE 18714 presented at the 
1989 SPE/IADC DriUing Conference, New Orleans, LA Feb.-March 1989. All of the 
above papers are incorporated herein by reference. Each technique can be readily 
modified by one of ordinary skill in the art, with the benefit of this disclosure, to 
utilize the NMR-derived bound water porosity in shales to provide an estimate of pore 
pressure in the shale and in a more permeable zone below the shale. 

[0018] For example, resistivity is often utilized to calculate water porosity. In 
normally pressured shales (hydro-pressure regime), there is a decrease of conductivity 
with depth of burial, due to compaction and associated water loss, that is, loss of 
presumably bound water due to decreasing porosity. This decreasing of shale 
porosity forms a "normal compaction line" trend. Generally, such "normal line" can 
be established for a given area, depending on various factors this line may shift right 
or left, but the slope remains essentially the same. Shale-reading points that fall to the 
left (higher conductivity area) of the normal compaction line are zones of abnormally 
high pressure. Step-by-step procedure for applying the resistivity approach is 
described, for example, by Matthewes W.R., and Kelly, J. "How to predict formation 
pressure and fi^acture gradient," The Oil and Gas Joumal, February 1967, the contents 
of which is incorporated by reference. Utilizing NMR-based bound water in such a 
method is straightforward, and has an advantage over the prior art in that salinity 
variations, fi^actures, and the like which may give misleading information fi-om 
resistivity logs will not be apparent in the NMR-based porosity data. 

[0019] In normally pressured shales sonic (travel time) logs shows a decrease in 
interval transit time with depth of burial. Sonic logs are typically utilized to estimate 
porosity of a formation. Since interval transit time is a fimction of porosity, a 
decrease in time is associated with shale porosity decreasing with depth. When depth 
is plotted on a linear Y-axis versus shale interval transit time on the logarithmic X- 
axis, a straight line can be drawn trough the normal pressure points. Points that fall to 
the left of the normal line are zone of abnormally high formation pressure. See SPE 
paper 62981, October 2000, and SPWLA paper BBB, June 2001, which are 
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incorporated herein by reference. This approach is less dependent on water 
conductivity, but strongly influenced by natural and drilling-induced fractures that 
may drive calculated porosity values too high. 

[0020] In a preferred embodiment, a method for while-drilling measurement of 
the formation pressure is based on bound water determinations from motion 
insensitive Ti (longitudinal polarization time) recovery process as the primary source 
of spectral porosity data. In another embodiment, the spectral porosity is determined 
from the T2 logging- while drilling measurement. In yet another embodiment, the Ti 
data from wireline logging tool is used to estimate formation pressure. In still another 
embodiment, T2 data from the wireline NMR logging tool is utilized to determine 
spectral porosity. In yet another embodiment, spectral porosity from both insensitive 
Ti and from the T2 data is collected, wherein Ti data is obtained from logging-while 
drilling measurements and the T2 data is collected either by logging-while drilling 
measurements or by wireline logging or by both in combination. In yet another 
embodiment, spectral porosity from both insensitive Ti data are collected, wherein Ti 
data is obtained from logging- while drilling measurements, and this data is combined 
with one or more of T2 data, sonic log data, and/or resistivity log data, all of which are 
beneficially collected by logging-while drilling. The data is interpreted to deduce 
increases in the water content of shales, which increases are indicative of abnormal 
pressure. In preferred embodiments of the above described methods, the data is 
automatically interpreted and the driller is advised of the probability of abnormal 
pressure substantially in real time, as the NMR logging tool enters the shale. Since 
NMR is a substantially lithology-independent measurement, the presence of various 
minerals associated with shale diagenesis or change in facies that adversely affect 
prior art methods of determining shale porosity will effect neither Ti nor T2 data 
obtained in accordance with the present invention. 

[0021] The method(s) used in accordance with this invention may be used to 
enable: (1) detection of the top of overpressure formations at the early stages of 
formation penetration, and (2) application of necessary changes to drilling program 
and casing design in real-time. In particular, unlike any other known logging method 
(e.g., resistivity, acoustic, density), the real-time Ti measurements enable handling of 
complex reservoir conditions including salt proximity, fracturing associated with 
faulting and mixed rock lithology. Notably, the disclosed method(s) and system are 
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extremely valuable in the areas where a shallow gas presence may cause significant 
drilling accidents (e.g., Gulf-of-Mexico, North Sea, Caspian Depression basins etc.). 



BRIEF DESCRIPTION OF THE DRAWINGS 

[0022] This invention may best be understood with reference to the attached 
drawings, in which: 

Fig. 1 is a plot of normal compaction trend cind resistivity LWD data versus 
depth for shale intervals. 

Fig. 2 is a plot of normal compaction trend and acoustic LWD data versus 
depth for shale intervals. 

Fig. 3 is a comparison of resistivity LWD data and acoustic LWD data from 
the same well. 

Figs. 4A-4D illustrate the impact of drilling on the NMR measurements. In 
particular, Fig. 4A is a schematic cross-section diagram of a NMR tool in a borehole 
showing the sensitive region as a circle around the tool within the borehole. Fig. 4B 
illustrates the effect of lateral tool movement on the sensitive region. Fig. 4C shows 
magnetization as result of a saturation pulse. Fig. 4D shows the effect of lateral tool 
motion on the readout sequence following the saturation pulse. 

Fig. 5 illustrates an example pulse sequence in accordance with a preferred 
embodiment. 

Fig. 6 represents a NMR- WD Ti log from a deep-water exploration well. 

Fig. 7 represents a NMR- WD Evaluation log from the same reservoir as 
shown in Fig. 6. 

Fig. 8 illustrates comparison between NMR- WD porosity data and acoustic 
log data. 

Fig. 9 illustrates a method for pressure prediction using Equivalent Depth 
Method. 
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Fig. 10 illustrates diametrical cross-sectional views of the NMR logging 
apparatus. 

Fig. 1 1 presents published data on T2cut-ofr values to separate CBW and 
effective porosity. 

Fig. 12 illustrates an exemplary T2 porosity distribution. 

Fig. 13 illustrates NMR response in a rock sample. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0023] Section A below describes the preferred embodiments of the method(s) 
used in accordance with this invention for the NMR-based porosity measurement and 
subsequent abnormal formation pressure determination. Section B discusses pore 
pressure determination from NMR data. Section C describes an example embodiment 
of a NMR LWD tool suitable for in situ abnormal pressure determination in 
accordance with the invention. For a comprehensive review of the NMR technology, 
including logging methods and various tool designs, the interested reader is directed, 
for example, to the book by Coates et al. entitled "NMR Logging: Principles and 
Applications" distributed by Gulf Publishing Company, and incorporated in its 
entirety by reference for background. 

(A) Porosity From NMR data 

[0024] It is well known in the field of NMR measurements of earth formations 
surrounding a borehole to use a downhole static magnetic field Bo to align the 
magnetic moment of spinning hydrogen (H^"^) protons in the formation in a first 
direction, the direction of the static magnetic field Bq. In order to establish 
equilibrium, the hydrogen protons must be exposed to the polarizing field for a 
multiple of the characteristic relaxation time Tj. Then, the magnetic component of a 
radio frequency (RF) electromagnetic wave pulse, which is polarized in a second 
direction orthogonal to the static field Bo, is used to tip the protons to align them in a 
direction that is orthogonal to both the first and the second direction. This initial RF 
pulse is thus called a 90"^ pulse. 

[0025] Following the 90° pulse, the protons in the formation begin to precess 
about the axis of the first direction. As a result, the protons produce an oscillating 
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magnetic field. However, due to inhomogeneities in the static magnetic field and 
irreversible molecular processes, the protons begin to de-phase, which causes the 
induced signal to decay. Nevertheless, the dephasing due to inhomogeneities in the 
static magnetic field is partially reversible. Therefore, by applying a 180*^ RF pulse, 
the instantaneous phases are reversed such that the protons gradually come back into 
phase, thus rebuilding the induced signal. After the signal peaks at the time when the 
protons are back in phase, the signal will then begin to decay due to dephasing in the 
opposite direction. Thus, another 180° RF pulse is needed to again reverse the 
instantaneous phases and thereby rebuild the signal. By repeating a series of 1 80"* RF 
pulses, the signal is periodically rebuilt after each dephasing, although each rebuilding 
is to a slightly lesser peak amplitude due to the irreversible molecular processes. This 
sequence is known in the art as the Carr-Purcell-Meiboom-Gill (CPMG) pulse 
sequence. Each rebuilding of the signal in this manner is called a spin echo, and the 
time constant associated with the decay of the spin echo amplitudes is called the 
transverse relaxation time, T2. Eventually, the irreversible processes prevail, such that 
no fiirther rephasing is possible and the signal dies out completely. 

[0026] In wireline logging, NMR-derived porosity is traditionally based on T2 
measurements. Notably, it was shown by Kleinberg et al. in "Nuclear Magnetic 
Resonance of Rocks: Ti vs. T2," paper SPE 26470, presented at the 68th Annual 
Technical Conference and Exhibition of the Society of Petroleum Engineers, Houston, 
TX, Oct. 3-6, 1993, that relaxation time measurements entirely based on Ti are not 
suitable for a moving wireline tool. But in accordance with the present invention this 
conventional understanding is challenged, and significant benefit is derived by using 
NMR relaxation time measurements determining longitudinal relaxation times Ti, 
instead-of transversal relaxation times T2 that are typically used by a wireline tool. In 
particular, the method(s) of this application take advantage of the magnetic field 
gradient, which can be approximated in the proximity of the sensitive volume as a 
linear fall-off of the magnetic field strength (and also of NMR resonance fi^equency) 
in the radial direction. Compared to T2 measurements, Ti acquisitions are more time 
consuming and limit logging speeds to 1-3 ft/min. For LWD, however, this speed 
limit agrees well with the maximvun rates of penetration (ROPs) to be expected (also 
1-3 ft/min.). The inherent robustness of a Ti measurement makes it a preferred 
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vehicle in accordance with this application for a while-drilling measurement of 
porosity and free-fluid volume. 

[0027] One problem with T2 measurements is that they are adversely effected by 
typical lateral and vertical tool motions that occur during logging-while-drilling 
operations. Figs. 4A-4B illustrate the impact of drilling on a regular data acquisition: 
as seen from above, the sensitive volume of the tool randomly moves with the tool 
motion. Fig. 4A shows a schematic cross section of the NMR tool in a borehole. The 
sensitive region is a cylinder around the tool within the formation. A T2 
pulse/measurement sequence started at any given point in time does not stay focused 
on a specific measurement volvime. Fig. 4b shows how random lateral tool movement 
displaces the sensitive region within the duration of a regular CPMG pulse sequence 
(100-500ms). Consequently, the returned echo amplitudes appear to decay much 
faster compared to a static measurement. A simple bound fluid vs. free fluid 
determination requires at least 100 ms measurement time. Given average lateral 
velocity and displacement conditions, the combined affects of dephasing and loss of 
measurement volume act to suppress the echo amplitudes within a few tens of 
milliseconds. 

[0028] Unlike T2 measurements, which are sensitive to magnetic field changes 
and the associated dephasing effects, the Ti recovery process is immime to small field 
fluctuations caused by lateral motion of the tool during drilling. Figs. 4c-4d illustrate 
Ti relaxation model. As shown in Fig. 4c, the Ti relaxation mechanism is 
characterized by saturation pulse that zeros nuclear magnetization within a large 
volume (e.g. 30cm in radius) shown as a gray band in Figs. 4c and 4d. Then, after a 
variable recovery time period that allows polarization of protons in the volume of the 
formation exposed to the magnetic field, the build-up of magnetization is monitored 
by fast readout sequences. Since the saturated volume of the formation is large 
compared to the lateral motion of the tool, the readout sequence is unaffected by the 
lateral tool motion. Accordingly, in a preferred embodiment of the invention, Ti 
relaxation spectra is used to derive formation porosity because Ti reflects the intrinsic 
relaxation mechanisms with a minimum of instrumentation artifacts. 

[0029] Finally, in accordance with the present invention the produced Ti 
relaxation curve is used to derive petrophysical properties of the formation or shale, as 
known in the art. In particular, the resultant T| relaxation curve is processed to 
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extract the dominant Ti relaxation modes, from which amounts of bound water, free 
water and hydrocarbons are estimated. The characteristic Ti times of the surface- 
wetting phase can also be used to estimate formation pore size distributions and 
formation permeability, especially in sandstone. 

[0030] The petrophysical model for Ti -based porosity is described next. The 

M 

NMR-measured total porosity is defined as ^f^^j^ = — , where Afioo%(0) is the 

^1000/. (0) 

measured magnetization for 100% bulk water with the same sensitive volume and Mq 
is the macroscopic net magnetization. The following equation defines a relationship 
between the longitudinal relaxation time Ti and the magnitude of magnetization Mz 
along the z axis, when protons are exposed to the Bo field for time Tw: 

-Tw 

M^iTJ = M,(l^e'^ ). 

[0031] To avoid underestimating Mo, fiiU polarization or nearly fixll polarization 
(95%) is desired. To achieve 95% polarization, Tw should be equal to three times Ti . 
It should be noted that the longitudinal relaxation time Ti has two components given 

by the expression: ^ = + , where Tm and Tis are the corresponding B - bulk 

and S - surface components, hi shales that are used for abnormal pressure 
characterization only surface component are typically present, and Ti could be 

expressed as:^ = -^-=pj^j , where p is the surface relaxivity and SA^ is 

^1 ^15 J PORE 

surface-to-volume ratio of pore system. 

[0032] Two data acquisition modes are developed using the petrophysical model 
for the Ti -based porosity. The first mode is called Ti -Logging (Ti); it is best suited 
for while-drilling operations because of low sensitivity to lateral and vertical tool 
motion. The second mode used in specific embodiments is T2-Logging (T2). It 
closely matches wireline logging (e.g., MRIL-Prime™ tool, available from 
Halliburton Energy Services, Inc., Houston, Texas) and thus is more suitable for 
sliding, wiping and tripping measurements. In accordance with a specific 
embodiment, either a logging-while-drilling or wireline tool may be used to acquire 
NMR data using these two modes. Preferably, the logging tool will alternately switch 
between the Ti and T2 modes. The transition between the data acquisition modes may 
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be triggered, for example, by elapsed time, by counting measurements, or can be tied 
to the outputs of accelerometers/magnetometers that distinguish drilling and non- 
drilling conditions. For instance, the while-drilling tool may be programmed to 
continuously toggle back and forth between the Ti and T2 modes. In postprocessing, 
drilling and non-drilling periods can be identified and the invalid T2 data recorded 
during drilling can be discarded. 

[0033] Table 1 summarizes the data acquisition parameters for the wireline tool 
(using the MRIL-Prime™ as base comparison), for the Ti mode (while-drilling 
logging) and T2 mode (after-drilling logging). T2 parameters closely resemble those 
of a wireline tool, except for slightly lower fi-equencies, field gradients and lower 
logging speeds. T2 speeds should be preferably limited to about 300 ft/hr if 
hydrocarbon typing (i.e., gas detection) is a logging objective. 



Table 1 - Comparison of Tool Parameters for 
MRIL-PrimeTM, MRI-LWD in Tj mode, and MRI-LWD in T2 mode 


Tool parameters 


MRIL-PrimeTM 
wireline 


MRI-LWD 
Ti mode 


MRI-LWD 
T2 mode 


Frequency Range 


580-750kHz 


500kHz 


580-750kHz 


Frequency Spacing (min) 


12kHz 


N/A 


12kHz 


No. of Frequencies (max) 


9 


1 


2 


Field Gradient 


17-32G/cm 


14G/cm 


~14G/cm 


Echo Spacing (min) 


0.5ms 


0.5ms 


0.5ms 


Measurement Mode 


T2 


T, 


Ti and T2 



[0034] The Ti log (Ti) used in a preferred embodiment is robust against typical 
lateral and vertical tool displacement due to its reliance on Ti measurements. Fig. 5 
illustrates the preferred mode for data acquisition during drilling based on a Ti 
recovery mechanism. The pulse diagram in Fig. 5 starts with a saturation pulse. This 
pulse reaches a large volume (gray band in Figs. 4c and 4d). The nuclear 
magnetization is instantly zeroed within this band during the pulse, and then allowed 
to recover to its equilibrium condition. After a timed recovery delay, a readout 
sequence consisting of two pulses and a single echo acquisition converts the 
recovered magnetization into an amplitude-modulated signal. The readout takes only 
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0.5 ms, which is short enough to suppress dephasing effects. Equally important, the 
signal beneficially comes from a thin-walled cylinder (about 1 mm thick). As long as 
this cylinder is contained in the much larger volume reached by the saturation pulse, 
the measurement is valid. The sequence is preferably repeated with variable recovery 
times, i.e., 1, 3, 10, 30, 100, 300, 1000 and 3000 ms. This allows various treatments 
of the data collected, which can provide different information as needed. 

[0035] In another embodiment, a method of obtaining Ti measurements comprises 
the steps of generating at least one radio frequency pulse covering a relatively wide 
range of frequencies to saturate the nuclear magnetization in a cylindrical volume 
around the tool; transmitting a readout pulse at a frequency near the center of the 
range of covered frequencies, the readout pulse following a predetermined wait time; 
applying at least one refocusing pulse following the readout pulse; receiving at least 
one NMR echo corresponding to the readout pulse; repeating the above steps for a 
different wait time to produce a plurality of data points on a Ti relaxation curve; and 
processing the produced Ti relaxation curve to derive petrophysical properties of the 
formation. 

[0036] In another embodiment, NMR data can be obtained by, for example, the 
method described in United States Patent No. 6,232,778, the contents of which are 
incorporated herein by reference. NMR measurements, in general, are accomplished 
by causing the magnetic moments of nuclei in a formation to precess about an axis. 
The axis about which the hydrogen nuclei precess may be established by applying a 
strong, polarizing, static magnetic field to the formation, such as through the use of 
permanent magnets. This field causes the proton spins to align in a direction parallel 
to the applied field (this step, which is sometimes referred to as the creation of 
longitudinal magnetization, results in the nuclei being "polarized"). Polarization does 
not occur immediately, but instead grows in accordance with a time constant Ti, and 
may take as long as several seconds to occur. After sufficient time, a thermal 
equilibrium polarization has been established. The 6,232,778 patent discloses 
methods for determining the bound fluid volume of a formation utilizing nuclear 
magnetic resonance (NMR) techniques in which less than fiiU polarization occurs and 
in which less than a complete NMR distribution is acquired. The effect of the 
polarizing static magnetic field is shortened in time by applying a shortened wait time 
between NMR measurements so that only the bound fluid of the formation is 
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polarized. The shortened wait time is effectuated by early application of the 
oscillating magnetic field to the formation which, in conjunction with a limited 
number of refocusing pulses, induces signals in the formation that are measured by 
the NMR tool. The peak amplitudes of these signals corresponds to the clay-bound 
and free-fluid volumes. 

[0037] In accordance with a preferred embodiment, the primary data record in Ti 
mode is a series of echo amplitudes associated with increasing recovery times and 
therefore of monotonically increasing amplitude. It is possible to collect more than 
one echo per readout sequence, but typically the first echo suffices. The readout 
sequence is phase-cycled similar to a CPMG (Carr-Purcell-Meiboom-Gill) sequence 
and the resultant echoes are co-added to eliminate instrumentation artifacts. As usual, 
the data points are corrected for environmental factors such as ambient temperature. 
The data points are fitted to a multi-exponential Ti model, as known in the art. The 
Ti times to be extracted ("bins") preferably are spaced on a logarithmic grid, such as 
the following example sequence: 1, 3, 10, 30, 100, 300, 1000, and 3000 ms. Although 
these happen to be the saturation-recovery times as well, such a one-to-one 
relationship is not necessary. Finally, in accordance with the method, the best-fit 
amplitudes for each bin are computed. 

[0038] Next, the porosity {^^^ra ) is computed as the sum of all amplitudes for all 
Ti bins. The resulted porosity is total porosity of water, oil and gas. In a specific 
embodiment, the Free Fluid Index (FFI) may then be computed as the sum of the 
amplitudes for Ti greater than 50 ms (bins 100 ms and higher). Then, a difference 
between the total porosity (p^j^j^ and the FFI is taken to determine the bound water 
porosity. In non-consolidated formations, which are characteristic to most drilling 
environments, only water-saturated shales are present, thereby the above-difference 
represents porosity associated with clay bound water (CBW). 

[0039] In accordance with another embodiment, the primary data record in T2 
mode is CPMG spin-echo amplitudes equivalent to those collected with the wireline 
tool using multiple wait times. The frequencies are advantageously, but not 
necessarily, programmable around a center frequency of about 500kHz, with a 
minimum echo spacing (Te) of 0.5 ms. The echo trains are phase-cycled and fitted to 
a T2 model of decaying exponentials. Multiple wait times allow computation of echo 
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differences, which are analyzed for the presence of oil and gas signatures. The 
procedure has been described in detail in the wireline logging literature and need not 
to be repeated in detail here. Final outputs include, for example: (1) T2 Porosity, 
corrected for gas hydrogen density where appropriate; (2) T2 Free Fluid, also 
corrected for gas hydrogen density; (3) T2 Capillary Bound Water, which is derived 
from the T2 distribution; (4) T2 Clay-Bound Water, which is also extracted from the 
T2 distribution. Fig. 12 illustrates a distribution of these porosity components on a 
typical T2 curve, with break points identifying CB W, BVI and FFI porosities. 

[0040] With reference to Fig. 12, several methods are known in the art to 
determine the break point on a T2 spectral porosity distribution curve. In one 
embodiment, a spectral BVI (SBVI) method may be used to determine the break 
points in the T2 spectral porosity distribution. The method is based on the recognition 
that each pore can contain both free and bound fluids More specifically, each pore 
size seen in the 100% brine-saturated spectra is assumed to contain some bound 
water. This is defined by a weighting ftmction W(T2,i), where 0 ^(T2,i) ^, that 
defines the fraction of bound water associated with each pore size. Various methods 
for obtaining the weighted functions have been described in the wireline logging 
literature and need not to be repeated in detail here. For more information on the 
spectral BVI see, for example, Coates, G., Xiao, Lizhi, and Prammer, M.G., NMR 
logging principles and applications; Halliburton Energy Services, 1999; 

[0041] In another embodiment of the present invention, cutoff BVI (CBVI) 
method is used to determine the break points on the T2 spectral porosity distribution. 
The method is based on a fixed T2 value (T2cutoff) that divides the T2 distribution into 
two components, one consisting of pore sizes containing bound water and the other 
consisting of pore sizes containing free fluids. In the absence of laboratory data, 
default values based on lithology are used for T2cutoff. A T2cutoff of 33 ms is used for 
sandstones, and 90-100 ms for carbonates; these values, however, are also affected by 
several other factors, such as pore- wall chemistry, minor paramagnetic or 
ferromagnetic components, texture, pore throat to pore body ratios, and other factors 
not well understood. In the same manner, applying a fixed T2 cut-off value of 3 ms 
identifies a micro-porosity associated with clay-sized particles. As seen from Fig. 13, 
which shows a sample of rock from Matagorda Island, the T2cutofr values of 3ms and 
33ms separate clay bound and capillary/free volumes correctly. 
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[0042] In yet another embodiment, Tacutoff can be detemiined in the laboratory 
with NMR measurements on core samples. Core samples are analyzed for NMR 
characteristics at two saturation conditions, Sw = 100% and (after establishing the 
appropriate value of saturation from a capillary-pressure curve, or directly 
desaturating the sample to the appropriate capillary pressxire) Sw = irreducible. 
However, this approach provides measurements of total amoxmt of irreducible fluids 
and does distinguish between clay- and capillary bound components. Various clay 
minerals exhibit wide range of relaxation times, some exceeding the commonly used 
3 ms cut-off For example, Prammer M.G., et al., in "Measurements of Clay-Bound 
Water and Total Porosity by Magnetic Resonance" the Log Analyst, 1996 37,#6,61 
reported T2 models at 2, 5 and 16 ms, respectively for clay-bound water in illite, 
chlorite and kaolinite, other published material gives different values. 

[0043] It should be further noted that when certain parts of T2 spectra are 
associated with clay or with intergranular porosity it is often assumed that pore sizes 
and distribution are proportionally related to the size and distribution of particles 
composing this rock. Thus, in clastic 100% water saturated environments good 
"clean" peak at large T2 values is associated with good well-sorted pores, which is 
typically representative of a clean well-sorted medium-to-coarse grained sandstone. 
When, however, observed peak is shifted towards the faster T2 components and 
become more widely "sagged", there is a good ground for interpreting such rock as a 
shaly fine-grained poor sorted sandstone. If that is true, then a one fixed cut-off of 3 
ms separating our microporosity associated with clay-sized particles might not be 
valid; especially in the cases represented by a uni-modal T2 distribution. Accordingly, 
in a preferred embodiment, minimum T2 value is used to separate CBW and free-fluid 
(or effective) porosity for bi-modal distribution. 

[0044] Furthermore, since NMR interpretation always utilize a various T2cut-off 
methods to discriminate free- and capillary-bound fluids, it is natural to investigate 
whether an additional cut-off exists which will separate clay bound water. Over past 
years several lab studies has been undertaken in that direction; some published results 
are summarized in the Fig. 1 1 . These result show that in general clay minerals exhibit 
the decrease in their relaxation times (T2) from kaolinite to chlorite, illite and 
smectite. The surface relaxivity of a clay is somewhat a fiinction of clay's iron 
content, and hence mineralogy. Kaolinite would be expected to mainly affect the 
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surface-to-volume ratio, whereas chlorite due to its iron content would also affect the 
surface relaxivity. The position of T2cut-off associated with the clays also depends on 
their distribution, i.e. rock diagenesis. But all exact values are descriptive rather than 
diagnostic, and there is no way that one single T2cut-off value will adequately describe 
the separation of clay bound porosity and free fluid porosity. 

[0045] Fig. 6 displays a small portion of a Ti Log (Ti), in accordance with this 
invention, where all data points were acquired during the actual drilling operation of a 
deepwater exploration well. From left to right, the tracks contain the following 
curves: Track 1 shows LWD natural gamma ray; Track 2 displays LWD resistivity 
logs for four different depths of investigation (i.e., extra shallow, shallow, medium 
and deep); Track 3 shows sonic, neutron and density porosities; and Track 4 presents 
NMR-WD porosity (the free fluid portion is shaded in black, and the bound fluid 
portion is shaded in light gray). From the NMR log it is observed that the reservoir is 
clearly identifiable by the increase in free fluid volume. This reservoir contains gas, 
which explains the drop in neutron porosity and the very high reading in sonic and 
density porosities. The ^;va//m porosity, however, undercalls true porosity due to the 
depressed hydrogen density of methane gas. 

[0046] For comparison. Fig. 7 shows an Evaluation Log acquired while wiping 
over the reservoir zone. Again, Track 1 shows natural gamma ray; Track 2 presents 
the reference MRIL-Prime™ wireline log (total porosity, capillary-bound fluid and 
clay-bound fluid); Track 3 shows T2 effective porosity, T2 free fluid and T2 capillary- 
bound fluid. For comparison, Ti porosity acquired while drilling is overlaid as dotted 
line, hi Track 4, T2 effective porosity and free fluid are shown corrected for the 
hydrogen index of methane gas for these reservoir conditions. As seen, the reservoir 
all results are very consistent. 

[0047] The agreement of total porosity obtained during Ti and effective porosity 
obtained during T2 indicates two facts: (a) the reservoir is clean sand with little clay 
contents, and (b) the invasion profile did not substantially change between drilling, 
wiping, and wireline logging. In particular, above, but also below the reservoir zone, 
there is total porosity, but no or little effective porosity. Li these intervals, fast 
relaxing clays with T2 less than 1 ms are prevalent. Such fast relaxation times are not 
resolved by the L2 ms echo train used for the wiping log. The LWD measurement 
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(echo spacing 0.5 ms), however, can detect fast clay components resulting in a correct 
total porosity reading. 

[0048] Another known method for calculating bound water using NMR is the 
fixed cutoff method, which is described in U.S. Pat. No. 5,389,877, the contents of 
which are incorporated by reference. That reference describes using NMR 
techniques, in which a short train of spin echoes (i.e., j echoes) are used to obtain a 
sharp cutoff (i.e., the fixed cutoff relaxation time) that can be used to determine the 
producible volume in a borehole by measuring bound fluid and subtracting it fi-om 
total porosity. They utilize a significantly shortened wait time tr, where tr is the time 
fi-om the last 180° pulse of one NMR experiment until the application of the 90° pulse 
to start the next NMR experiment— where only single pulse sequences are utilized 
(e.g., while the drill pipe is moving), the shortened wait time is defined from the time 
the spins enter the static magnetic field. The shortened wait time is selected so that 
the static magnetic field is only appUed for a short period of time before the next 
NMR experiment begins. This results in polarizing only the bound fluid in the 
formation, rather than the bound and free fluid, so that the measured amplitude 
represents the actual bound fluid volume of the part of the formation being analyzed. 

(B) Pore Pressure from NMR data 

10049] The estimation or determination of pore fluid pressure is a major concern 
in any drilling operation. The pressure applied by the column of drilling fluid must be 
great enough to resist the pore fluid pressure in order to minimize the chances of a 
well blowout. Yet, in order to assure rapid formation penetration at an optimum 
drilling rate, the pressure applied by the drilling fluid column must not greatly exceed 
the pore fluid pressure. Likewise, the determination of horizontal and vertical 
effective stresses is important in designing casing programs and determining pressures 
due to drilling fluid at which an earth formation is likely to fracture. 

10050] The commonly-used techniques for making pore pressure determinations 
have relied on the use of overlay charts to empirically match well log data to driUing 
fluid weights used in a particular geological province. These prior art techniques are 
semi-quantitative, subjective and unreliable from well to well; none is soundly based 
upon physical principles. 
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[0051] Effective vertical stress and lithology are the principal factors controlling 
porosity changes in compacting sedimentary basins. Sandstones, shales, limestones, 
and others, compact at different rates under the same effective stress. An effective 
vertical stress log can be calculated from observed or calculated porosity for each 
lithology with respect to a reference curve for that lithology. Calculating pressure, 
and detecting abnormal pressure, is particularly beneficial while drilling through 
shales. Drilling through overpressured shales will not cause a large volume of low 
density fluid to enter the wellbore and displace drilling mud, because the permeability 
of the shales is too low. However, if a high permeability region with abnormal 
pressure, for example a typical sandstone reservoir, is encountered with 
underweighted drilling mud, low density fluid will enter the wellbore at increasing 
rates and the well might become uncontrollable. 

[0052] Therefore, once the porosity and bound water volume of the shale 
formation has been determined, various techniques may be used to determine pressure 
deviations from the normal compaction trend. Such known techniques are described 
by (1) Hottmann et al. in " Estimation of Formation Pressures From Log-Derived 
Shale Properties" JPT {Journal Petroleum Technology), June 1965; (2) Eaton in "The 
Equation for Geopressure Prediction From Well Logs," paper SPE 5544, 1975 Annual 
Technical Conference, Dallas, Sept. 28 - October 1 ; (3) Rasmus et al. " A Framework 
To Estimate Pore Pressures in Real-Time," February 1990. Measurement- While- 
Drilling Symposium, Louisiana State University, Baton Rouge, LA; (4) Holbrook et 
al. in "A Petrophysical/Mechanical Mathematical Model for Real-Time Wellsite Pore- 
Pressure/Fracture-Gradient Prediction, Paper SPE 16666 presented at 1987 SPE 
Annual Technical Conference , Dallas, September; and (5) Bryant in "A Dual Shale 
Pore Pressure Detection Technique" paper SPE 18714 presented at the 1989 
SPE/IADC Drilling Conference, New Orleans, LA Feb.-March. All of the above 
papers incorporated herein by reference and are familiar to those of ordinary skill in 
the art. 

Example 1: Ben Eaton*s Method 

[0053] Ben Eaton's method for pore and fracture pressure analysis is one of the 
most widely industry-accepted techniques. Particularly, it has proven to be extremely 
valuable in the Gulf-of-Mexico area. Two key equations, developed by Eaton, after 
modifications are as follows: 
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p =p - 



-0.465*1^^^' 
V Vnmra j 



E 



where: Pp - pore pressure or formation pressure gradient in psi/f 
Po - overburden pressure in psi/f; 

0.465 - water gradient, suggested by B. Eaton for Gulf-of-Mexico; 

^MRN - normalized NMR-porosity for given depth; 

(f^MRA - apparent (measured) NMR-porosity for given depth; 

E • Eaton coefficient, which was suggested to be 1 .2 when resistivity data used 
and 3.0 for acoustic data; however later research indicate that depending on reservoir 
conditions, lithology, and the like this coefficient may vary fi-om 1.2 to 5.0. 
Advantageously, an Eaton coefficient is determined for a region and utilizing NRM 
bound water porosity values. 

[0054] Once the pore pressure heis been computed the fi-acture pressure may be 



calculated as follows based on: = " — — + Pp 



^ where v- Poisson's ratio. 
\-v ) 

[0055] It should be noted that the overburden pressure at given depth {D[ft]) 

D 

equals cumulative weight of the rocks above: P^ = 0.4335 jpf,(D)dD where 0.4335 - 

0 

is a conversion coefficient of g/cc into psi/f; and = ^f^^iuPw + (1 - (I>nmra)Pm i-e. 
fimction of porosity; w - water and ni-matrix density. If no density information is 
available, as for example in the research conducted in Gulf-of-Mexico, studies allow 
calibrating changes in shale density vs. depth which could be approximated with the 
following exponents: pb = 0.2424*ln(D) + 0.23 for consolidated formations, and pb = 
0.1661*ln(D) H- 0.84 for unconsolidated formation. 

[0056] The overpressurized zones may be determined very accurately using Ti- 
derived porosity and Eaton's pressure analysis method. Fig. 8 presents comparison of 
acoustic LWD log for shale intervals and CB W trend fi-om the NMR T2 measurements 
in accordance with the preferred embodiment of the invention described above in 
Gulf-of-Mexico Matagorda Island. In this study, the acoustic logs are believed to 
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have provided accurate results and there is a very good correlation between the 
acoustic logs and NMR measurements. 

Example 2: Equivalent Depth Method 

[0057] Another known method for formation pressure analysis that can be used 
according to this invention is the Equivalent Depth Method. With reference to Fig. 9, 
the technique generally requires several steps: First, porosity ^az/w is plotted versus 
depth and a normal trend line is established. Second, a vertical line that intercepts the 
normal trend of NMR porosity is drawn. Third, depths at point of interest "A" and 
point of interception "B" are recorded. Fourth, a formation pressure at point "A", 
= - 0.535i)2 , is calculated. 

[0058] It should be noted that many pore and fracture pressure models require the 
use of shale values for compaction trend. Therefore, in accordance with a specific 
embodiment of the invention, the following parameter may be used as a shale 

CBW 

indicator on spectral porosity data: for consolidated formations > Ti cut-off, 

BFI 

and for non-consolidated formations > Ti cut-off* The statistical studies 

CBW + BVI 

of overpressured formations in the Gulf-of-Mexico suggest a value of the above- 
parameter of 0.85, however, it may vary between various formations. 

[0059] Other methods or correlations are known and may be developed, which 
correlate the bound water content, and/or the porosity, of shale layers with the 
likelihood of overpressure, and these other methods or correlations are expected to be 
usefiil with this disclosed method. See, for example, U.S. Pat. No. 4,981,037, which 
is incorporated by reference. Furthermore, robust porosity measurements could be 
used in previously accepted pore pressure models or independently as an Equivalent 
Depth method. 

[0060] This method enables recognition of the presence of gas and is not sensitive 
to induced fractures and/or changes in formation water salinity. It is possible to 
utiHze data obtained from the Ti measurements to deduce the presence of 
hydrocarbon-based hydrogen, for example by utilizing a different part of the spectral 
data, as is known in the art. The porosity data can be corrected for the presence of 
hydrocarbons, or alternately the data can be collected so that the contributions of the 
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hydrogen associated with hydrocarbons is not a pertinent factor. Ahematively, the 
presence of hydrocarbons in a shale zone may, in some regions, itself be an indicator 
of overpressure in the underlying formations. Therefore, recognition of complex 
overpressured zones in Black Shales associated with liquid hydrocarbon generation 
also could be performed. 

[0061] In sum, a novel NMR-based technique has been disclosed that allows 
determination of the normal compaction trend and calculation of overpressure in 
transition zones using Ti real-time measurements. Most importantly, the NMR signal 
gives a "spectral" or incremental porosity, which allows the interpreter to separate the 
clay-bound water signal from the signal from capillary-bound/free fluids. The method 
of detecting overpressure can be based on the comparison of the clay-bound water 
porosity of shale zones, preferably shale zones that are separated vertically in a well 
by at least one permeable zone. Advantageously, this at least one permeable zone, 
i.e., permeability of at least one millidarcy, is normally pressured. The detection of 
overburden can be the result of anomalies in the normal clay-bound water in shale 
versus depth for a given well, or it can be a numerical calculation of pore pressure or 
other indicators that is/are based on correlations that have been shown to be reliable 
for a region. For example in a well drilled with a water-based mud, the presence of 
induced or tectonic fractures will "mislead" resistivity and acoustic logs. NMR on the 
other hand will discriminate between a signal from drilling mud in fractures (which 
will show-up in a BVI part of spectrum) and a micro-porosity associated with shales. 
For example, in wells drilled in the Bazenov formation (West Siberia, Russia), the 
presence of liquid hydrocarbons was a direct indicator of overpressure. Resistivity 
Logs could not have been used due to relative absence of water in the shales. 
Acoustic Logs were confusing due to extremely complex lithology (mixture of 
carbonate, siliceous material, dispersed organic matter and shales). Neutron logs read 
all hydrogen - even in the tar sections. In that case, NMR measurements were the only 
ones that allow evaluation of the formation for presence of liquid media and hence 
overpressure recognition in black shale formations. 

[0062] A preferred novel NMR-WD-based technique that allows determination of 
the normal compaction trend and calculation of overpressure in transition zones, i.e., 
shale zones, uses Ti real-time measurements to calculate the bound water fraction of 
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shale porosity, and then incorporates this data into one of several known correlations 
to detect and calculate abnormal pressure. 

(C) The NMR LWD Tool 

[0063] In a preferred embodiment of the invention, a NMR tool having the 
mechanical strength and measurement sensitivity to perform NMR measurements of 
an earth formation while drilling a borehole is utilized. Fig. 10 illustrates an NMR 
WD tool 40 used and operated in accordance with a preferred embodiment of the 
present invention. NMR tool 40 comprises a nonmagnetic metal drill collar 10 that 
encloses a tubular, permanent magnet 12, which surrounds a mud tube 34 through 
which drilling mud may be pumped during the drilling of a borehole 36 into the earth 
formation 138. The tool 40 further comprises a stabilizer section 66, antenna 14, 
magnet section 12, capacitor compartment 32 with tuning capacitors 42, 
accelerometers 24, antenna driver 28, signal processor 30 and data transmitter 150. 

[0064] Permanent magnet 12 has a magnetic field direction Bq that is substantially 
perpendicular to the axis of the borehole. Collar 10 of a non-magnetic material, or of 
a substantially non-magnetic material, may beneficially surround magnet 12. At least 
one antenna 14 is positioned a distance from the permanent magnet, particularly 
outside the collar. Advantageously, a soft magnetic material is positioned in a 
predetermined relationship with the collar and the magnet that helps to shape the 
magnetic field of the tool 40. Due to the non-magnetic collar 10, the tool 40 can 
withstand the extreme conditions in the borehole environment while the borehole is 
being drilled. Motion management apparatus and method are beneficially employed 
to identify time periods when the NMR measurements can be taken without the 
accuracy of the measurement being affected by the motion of the tool or its spatial 
orientation. 

[0065] While characteristics of a tool in a preferred embodiment have been 
presented, other tools are known and can be used to characterize formation properties 
in the process of predicting overpressure, including those tools described in U.S. 
Patent No. 5,923,167, U.S. Patent No. 6,11 1,408; U.S. Patent No, 6,173,793; U.S. 
Patent No. 6,204,663; U.S. Patent No. 6,242,912; U.S. Patent No. 6,242,913; U.S. 
Patent No. 6,255,819; U.S. Patent No. 6,268,726; U.S. Patent No. 6,316,940; U.S. 
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Patent No. 6,337,568; U.S. Patent No. 6,362,619, U.S. Patent No. 6,366,087; U.S. 
Patent No. 6,41 1,087; U.S. Patent No. 6,446,736; U.S. Patent No. 6,518,756; U.S. 
Patent No. 6,531,868; and U.S. Patent No. 6,541,969, which are incorporated herein 
by reference. Processing of NMR data is done in software, hardware or a suitable 
combination preferably capable of real-time data manipulation. Naturally, the 
selected processing mechanism can be programmed to implement the above described 
abnormal pressure zones automatically. In the altemative, porosity data from the tool 
can be displayed to a human operator who can make a decision regarding abnormal 
pressure zones based thereon. 

[0066] While the foregoing has described and illustrated aspects of various 
embodiments of the present invention, those skilled in the art will recognize that 
altemative components and techniques, and/or combinations and permutations of the 
described components and techniques, can be substituted for, or added to, the 
embodiments described herein. It is intended, therefore, that the present invention not 
be defined by the specific embodiments described herein, but rather by the appended 
claims, which are intended to be construed in accordance with the following well- 
settled principles of claim construction. 
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